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Abstract. Recent progress in constructing holographic models for QCD is discussed, concentrating on the bottom-up models 
which implement holographically the renormalization group flow of QCD. The dynamics of gluons can be modeled by using 
a string-inspired model termed improved holographic QCD, and flavor can be added by introducing space filling branes in 
this model. The flavor fully backreacts to the glue in the Veneziano limit, giving rise to a class of models which are called 
V-QCD. The phase diagrams and spectra of V-QCD are in good agreement with results for QCD obtained by other methods. 
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INTRODUCTION AND MOTIVATION 

Holographic QCD has been the topic of numerous stud¬ 
ies after the discovery of gauge/gravity duality. The stud¬ 
ied models can be roughly divided into two classes, the 
top-down and the bottom-up models. 

The top-down models are fixed constructions directly 
based on string theory and motivated by the original 
AdS/CFT duality. In particular, the dual field theory can 
usually be identified exactly, and the gravitational de¬ 
scription does not bring in any extra parameters: all pa¬ 
rameter of the gravitational action are uniquely fixed. 
The ultimate goal of this kind of constructions would 
be to to find an explicit holographic realization for the 
infrared (IR) physics of QCD, at large number of col¬ 
ors Nc- However, so far the available models have some 
shortcomings. Most importantly, when classical calcu¬ 
lations are reliable on the gravity side, the low-energy 
spectrum of the dual field theory has additional (Kaluza- 
Klein) states at the same mass scale as the states of QCD. 
Nevertheless, the observables computed in these mod¬ 
els are in good agreement with our knowledge of QCD 
(see [1, 2, 3] for concrete examples). 

In the bottom-up approach one usually takes only the 
main ideas from string theory. The holographic (often 
five dimensional) model is constructed “by hand” by 
picking the most important operators on the field the¬ 
ory side, and then writing down a natural action for 
their dual fields. This means that the coupling constants, 
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masses and/or potentials of the constructed action are not 
known, but need to be chosen by comparing the results to 
the field theory. Therefore the holographic description is 
necessarily effective and and a wide class of gauge the¬ 
ories is modeled instead of a definite single theory. Typ¬ 
ically the action is chosen such that it respects general 
symmetries and other features found in top-down con¬ 
structions. They are able to produce a very precise de¬ 
scription of QCD data (such as the meson spectrum) just 
by fitting a few parameters. Well-known examples of this 
type of models are [4, 5, 6]. Here we will discuss a spe¬ 
cific bottom-up approach which is more strictly based on 
string theory than is typical for bottom-up models. 

So far the majority of the literature on holographic 
QCD discusses flavor the’t Hooft or probe limit, where 
one takes Nc keeping Nf and the ’t Hooft coupling 
Sym^c fixed. Here we will, however, also consider the 
backreaction of the flavor to the glue, which is present if 
we take the Veneziano limit instead. It is defined by 

Nc-to° , X= fixed , 

Nf 

Nf^°°, x =— fixed. (1) 

Nc 

Backreaction means that the analysis will necessarily 
be technically involved. Taking it into account is, how¬ 
ever, important for several reasons. First, the number of 
light quarks V/ ~ 2... 3 is comparable to the number of 
colors Nc = 3 for ordinary QCD, so that the backreaction 
is expected to be strong. But there are also phenomena 
(e.g., phase transitions) which cannot be accessed in the 
probe limit at all. From arguments based on perturbation 
theory, it is clear that QCD has the so called conformal 



window: the theory has a nontrivial IR fixed point and 
chiral symmetry is intact for some Nf, Nc with not so 
small Nf/Nc (whereas in the ’t Hooft limit Nf/Nc —> 0). 

More precisely, in the Veneziano limit, the leading 
coefficient of the QCD beta function vanishes at x = 
11/2 = xbz- As X —7- xbz from below, the theory has an 
perturbative IR fixed point, i.e., a fixed point at a small 
^ (-^BZ — x) value of the coupling. But for small x the IR 
physics is expected to be similar to ordinary QCD, with 
confinement and chiral symmetry breaking. Therefore a 
phase transition must take place at some critical value Xc 
of the ratio x, which marks the lower edge of the confor¬ 
mal window. This transition has been called the confor¬ 
mal transition. It has been conjectured that right below 
the transition the IR fixed point is still visible as the the¬ 
ory displays “walking” (nearly conformal) behavior, i.e., 
the coupling constant stays approximately constant for a 
large range of energies. 

In the models discussed here, pure Yang-Mills is mod¬ 
eled by improved holographic QCD (IHQCD) [7]. This 
model is inspired by five dimensional string theory, but 
also contains a potential term which is not fixed a pri¬ 
ori. Flavor can be introduced in IHQCD by adding a pair 
of space-filling five dimensional D-branes [8]. Includ¬ 
ing full backreaction, by taking in the Veneziano limit, 
then leads to the V-QCD models [9]. Walking dynam¬ 
ics and the conformal transition have also been studied 
in closely related but simpler models called Dynamic 
AdS/QCD [10]. 

In these proceedings we will first discuss IHQCD and 
V-QCD in more detail. Then we will give an overview of 
the results obtained in the V-QCD models. 

IMPROVED HOLOGRAPHIC QCD 

Improved holographic QCD (IHQCD) [7] is a model 
for the Yang-Mills theory. It is loosely based on five 
dimensional noncritical string theory. It is defined in 
terms of the five dimensional action 

SmQCD=M^N'^ j d^x^/-del g 

( 2 ) 

where the first term defines the five dimensional Einstein 
gravity and 0 is the dilaton. 

We denote X = e‘^ and identify the ’t Hooft coupling 
as the vacuum solution X{r), with the vacuum being the 
extremum of the action (2). Here r denotes the fifth bulk 
coordinate. The Ansatz for the metric is 

ds^ = e^^^''\dx\^+dr^) (3) 

in the vacuum solution. The factor A is identified as 
the logarithm of the energy scale in field theory. In our 
conventions the ultraviolet (UV) boundary lies at r = 0, 


and the bulk coordinate runs from zero to infinity. The 
metric will be close to the AdS near the UV. 

The model deviates from the “top-down” derivation 
from five dimensional noncritical string theory as the 
dilaton potential Vg is fixed by hand rather than derived 
from string theory. Its asymptotic behavior is determined 
using as a guideline known properties of glue dynamics 
in Yang-Mills theory. In the IR (i.e., large values of the 
coupling X), potential which diverges as 

Vg^X^/^y^, (X^oo)^ (4) 

generates confinement, a mass gap, discrete spectrum 
and asymptotically linear glueball trajectories [7, 11]. 
Interestingly, this matches with the behavior predicted by 
string theory Vg ~ X^!^ up to the logarithmic correction, 
which is important to ensure exactly linear trajectories. 

In the UV (i.e., as A —> 0), the geometry is asymp¬ 
totically AdS and the coupling X has logarithmic de¬ 
pendence on r, similar the renormalization group (RG) 
flow of the coupling in Yang-Mills, when Vg is analytic 
at A = 0. The holographic beta function defined by the 
above dictionary, j3 = A'(r)/A'(r), then admits a Tay¬ 
lor expansion at A = 0 (with X evaluated on the vacuum 
solution). The coefficients in the Taylor expansion of Vg 
are chosen such that the holographic beta function agrees 
with the perturbative beta function of Yang-Mills up to 
two-loops. 

At intermediate couplings Vg must be fitted to Yang- 
Mills data. Most precise constraints are obtained by 
studying the thermodynamic functions [12]. Even a two 
parameter fit gives a very good description of the lattice 
data [13]. 

V-QCD 

A framework for adding flavor in holographic models has 
been introduced in [8]. It is realized through adding a pair 
of DA-DA branes which therefore fill the whole space in a 
five dimensional holographic model. The string stretch¬ 
ing between the branes gives rise to a tachyon field, 
which is dual to the qq operator on the field theory side, 
and therefore controls chiral symmetry breaking. The dy¬ 
namics of the tachyon T is governed by a tachyon-Dirac- 
Born-Infeld (TDBI) action, and taking an exponential 
potential, Vf{x) = exp(—T^), for the tachyon. The frame¬ 
work has been tested in the probe approximation [14] 
with encouraging results. 

The V-QCD models arise as the fusion of IHQCD with 
this framework for adding flavor, in the Veneziano limit 
where flavor fully backreacts to the glue [9]. The main 
degrees of freedom are therefore the “dilaton” X = 
and the tachyon T. The dictionary is as stated above: (the 
background value of) X is the’t Hooft coupling, and T is 
dual to qq. 



The terms in the V-QCD action which determine the 
vacuum solution are 


SyQCD=M^N^ J d^x^/-detg 

-M^NfNc J d^xVfiX,T) 

X ^-det {g^v + K{X)d^zdvT) . 


(5) 


The first term is simply the IHQCD action (2) and de¬ 
scribes the dynamics of the glue in the model. The sec¬ 
ond term is a generalized TDBI action, which describes 
the dynamics of the flavor. From the factors of Nc and 
Nf it is clear that both terms are leading in the Veneziano 
limit and there is indeed full backreaction. 

Notice that the form of the TDBI action is precisely 
known only in the probe limit. Therefore we have intro¬ 
duced a more general Ansatz in (5). For the tachyon po¬ 
tential we assume the form 

Vf{X,x)=Vfo{X)e-‘'^^^^" (6) 


where the dependence on the tachyon is the same as in 
the probe limit above. The potentials V/o(A), k{X), and 
a{X) would be power laws in the probe limit, but we have 
taken them to be arbitrary functions of A. 

In order to pin down the model, we need to specify 
the functions Vg{X), Vfo{X), k{X), and a{X). It turns 
out that the dilaton potential Vg{X) must obey the same 
asymptotic constraints which were found in IHQCD. For 
the other functions similar constraints are found, both in 
the IR and in the UV. 

In the IR (A —°°), the functions are constrained, 
among other things, by the asymptotics of the spectra at 
large excitation numbers, and by requiring that the ge¬ 
ometry has a “good” kind of IR singularity that uniquely 
fixes all IR boundary conditions. Consequently, we find 
the following asymptotics [15]: 

a(A) ~ const.; ff(A) ^ A^^/^\/logA , {X^°°). 

(7) 

Interestingly, these asymptotics again agree with the 
string theory prediction, up to the logarithmic correction 
in K'(A), in analogy to what was found for Vg{X) in (4). 

In the UV (A —^ 0) the potentials must approach con¬ 
stant values in order for the dimension of the qq opera¬ 
tor to be correctly reproduced, and for the geometry to 
be asymptotically AdS in the UV. In analogy to the ap¬ 
proach taken in IHQCD [7], we take the potentials to be 
analytic at A = 0 and match the coefficients of their ex¬ 
pansions with those of the perturbative beta function of 
QCD and the anomalous dimension of the quark mass 
7 = —d\ogmq/d\oggL. This matching is done in order to 
guarantee correct boundary conditions for the IR dynam¬ 
ics, which can be more reliably modeled by holography 
and is therefore more interesting in this model. 
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FIGURE 1. The phase diagram of V-QCD (at zero quark 
mass and temperature) as a function of x = Nf/Nc- 


The remaining freedom is the dependence of the po¬ 
tentials on A at intermediate values, which can be fitted 
to lattice and QCD data (neglecting the 1 /Nc corrections 
arising when one moves from the Veneziano limit to fi¬ 
nite Nc) in order to build a concrete model for QCD. Such 
a fit will be discussed in a future publication. The results 
below are based on potentials which obey the asymp¬ 
totics found above, but were not yet fitted to any data. 


PROPERTIES OF V-QCD 

In order to uncover the phase structure of V-QCD, we 
write an Ansatz where the fields A(r) and T(r) depend 
on the bulk coordinate r, and the metric has the form (3). 
The Ansatz is then inserted to the equations of motion de¬ 
rived from the action (5), which are consequently solved 
numerically. 

This procedure results in phase diagram (at zero tem¬ 
perature and quark mass) that is essentially universal 
and shown schematically in Fig. 1. The diagram has two 
phases separated by a phase transition at x = ~ 4: 

• In the region 0 < x < x^, the theory has chiral sym¬ 
metry breaking, and IR dynamics is similar to that 
of ordinary QCD. The corresponding background 
solution has nontrivial A (r), A{r) and T(r), with the 
tachyon diverging at the “good” IR singularity of 
the geometry. 

• In the conformal window, i.e., when Xe < x < 11/2, 
the theory flows to a nontrivial IR fixed point and 
is chirally symmetric. The background solution has 
zero tachyon T(r)=0 and nontrivial A(r) and A{r), 
giving rise to a geometry flowing to a nontrivial AdS 
fixed point in the IR [16]. 

As X — Xc from below, “walking” RG flow of the cou¬ 
pling constant is found. 

The location of the Banks-Zaks region near x = 11/2, 
and the existence of the perturbative fixed point is a 
consequence of matching the potentials to the perturba¬ 
tive RG flow of QCD in the UV. The conformal transi¬ 
tion and the walking regime, on the contrary, are gen¬ 
erated by the dynamics of the model in a nontrivial 
way [17, 18]. One can show that the conformal transition 
is related to a well-known instability of the tachyon field 
in the conformal window. This instability arises when the 
five dimensional bulk mass of the tachyon, evaluated at 
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FIGURE 2. Examples of phase diagrams of V-QCD at finite temperature on the (x, r)-plane (left) and at finite chemical potential, 
on the (T, /i)-plane with fixed x = 1 (right) [20, 21]. 


the IR fixed point, violates the so-called Breitenlohner- 
Freedman bound [19]. Transition triggered by the vio¬ 
lation of this bound quite in general leads to walking 
and a BKT/Miransky phase transition of infinite order. 
This transition involves a characteristic exponential scal¬ 
ing law. E.g., the chiral condensate vanishes as 

(gg) ^exp , (x^Xc^). (8) 

Phase structure at finite temperature and 
chemical potential 

The thermodynamics of the model can be studied as 
for IHQCD, by looking for solutions with a horizon, i.e., 
black holes. Rich structure of black holes, with one or 
two scalar hairs, was found in [20]. A typical phase dia¬ 
gram at zero quark mass on the (x, 7’)-plane is depicted 
in Fig. 2 (left). For low values of x, there is single first 
order deconfinement/chiral transition (solid blue curve). 
Near the walking regime, the chiral restoration happens 
at a separate second order transition (dashed blue curve). 
The crossover (dashed red curve) reflects the flow from 
the nearly conformal (walking) regime to the UV (run¬ 
ning) regime. The transition temperatures obey Miransky 
scaling as x —>■ Xc- 

Quark chemical potential can be added in the model 
following the standard recipe, i.e., by allowing for a 
nontrivial gauge field Ao(r) in the TDBI action in (5), and 
by setting ji — Ao{Q) [21]. An example of the resulting 
phase diagram is shown in Fig. 2 (right), where x = 
1. The red solid curve is the first order deconfinement 
transition, and the blue dashed curve is the second order 
chiral transition. As T —>^ 0 in the chirally symmetric 
phase, the IR geometry is asymptotically AdS 2 x and 
the entropy approaches a finite value instead of zero. 
These observations suggest that there is an instability. 


which may be linked to the color flavor locked or color 
superconducting phases of QCD. 


Spectra and the S-parameter 

The fluctuations around the vacuum solution of the V- 
QCD action have been studied at zero temperature [22, 
15] and at finite temperature [23]. Full analysis requires 
including flavor structure and gauge fields in the TDBI 
action (see [15] for details). There are both flavor singlet 
(corresponding to quark bilinear operators contracted 
with the unit matrix iNfxNf in flavor space) and non¬ 
singlet (corresponding to operators with the traceless 
Hermitean generators f“) fluctuations. The normalizable 
fluctuation modes are identified with meson states. 

The masses of the nonsinglet mesons are relatively 
easy to compute. We show the masses of the lowest mode 
in each sector (scalars, pseudoscalars, vectors, and axial) 
as a function of x in Fig. 3 (left). The dependence on x 
is weak for small x, but in the walking regime all masses 
go to zero obeying the Miransky scaling law. 

There is nontrivial mixing between mesons and glue- 
balls in the flavor singlet scalar sector. It has been con¬ 
jectured that the scalar singlet state contains a “dilaton” 
as X —> Xc, an anomalously light Goldstone mode due to 
the breaking of conformal symmetry. No such state was 
found in the analysis of the scalar singlet states of V- 
QCD. 

Another important quantity which can be determined 
from the fluctuations is the S-parameter. For a strongly 
interacting theory to be a viable candidate for the real¬ 
ization of technicolor, i.e., in order for the theory to pass 
the precision tests at FEP, this parameter must be much 
smaller than one. The common lore is that this parame¬ 
ter is suppressed in the walking regime (see, e.g., [24]), 
suggesting that walking technicolor theories are viable. 

The (normalized) S-parameter in V-QCD is shown 
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FIGURE 3. Left: The flavor singlet masses as a function of x in V-QCD. The mass of the lowest state at each sector is shown. 
Right: The S-parameter as a function of x in V-QCD at niq = 0 (the dashed red curve) and at niq = 10^® (the solid blue curve). 


in Fig. 3 (right). The blue solid (red dashed) curves 
have niq = 10^® (niq = 0), respectively [22, 15, 25]. 
We notice that the normalized S-parameter is enhanced 
rather than suppressed in the walking region (as x 
from below). It will be interesting to see if this behavior 
remains true after the potentials of the V-QCD action 
have been fitted to QCD data. Our results show that as 
a tiny quark mass is turned on in the conformal window, 
the S-parameter immediately jumps to a finite ^ (NfNc) 
value. Such a discontinuity is consistent with recent field 
theory analysis [26]. 

CONCLUSIONS AND OUTLOOK 

In these proceedings we have discussed the construc¬ 
tion of V-QCD and the most important results obtained 
with the model so far. This model takes into account 
the backreaction of the flavor degrees of freedom and 
therefore allows us to analyze the spectra as a function 
of X = Nf/Nc- The results are encouraging: the model 
mostly agrees, at qualitative level, with earlier literature 
on QCD. 

Several additional extensions and applications of the 
model are under study. These include 

• A study of the CP-odd sector, including the theta 
angle and axial anomaly of QCD. 

• An analysis of the dependence of meson masses and 
correlators on the quark mass. 

• Comparison of the thermodynamics to that of low 
temperature hadron gas. 

• Fitting the model to QCD data, for example Yang- 
Mills thermodynamics from the lattice, experimen¬ 
tal values of the QCD meson masses, and the equa¬ 
tion of state at hnite x from the lattice. 
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